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SUSY Signatures at LHC 

Prank E. Paige 

Brookhaven National Laboratory, Upton, NY 11973 USA 

Abstract 

The ATLAS and CMS Collaborations at the CERN Large Hadron Collider 
(LHC) have devoted considerable effort to the study of SUSY signatures and mea¬ 
surements. This talk provides an overview of what can be learned at the LHC if 
TeV-scale SUSY exists. 


1 Introduction 

SUSY is perhaps the most promising candidate for physics beyond the Standard Model; 
it also provides a good test of detector performance. The ATLAS [jl| and CMS p Col¬ 
laborations at the CERN Large Hadron Collider (LHC) have therefore devoted a lot of 
effort to studying SUSY signatures and measurements. This talk provides an overview of 
that work with emphasis on new results since the overviews in Ref. and 

The SUSY cross section at the LHC is dominated by the associated strong production 
of gluinos and squarks. If R parity is conserved, these decay into the lightest SUSY 
particle which escapes the detector, plus quarks, gluons, and perhaps other Standard 
Model particles. Thus, SUSY provides signatures containing at least jets and large missing 
transverse energy The LHC should be able to observe these signals for g and q masses 
up to about 2TeV with only 10fb“^ of integrated luminosity. 

The challenge at the LHC is not to discover TeV-scale SUSY (assuming it exists) but 
to make precision measurements of masses and other quantities. Since the decay prod¬ 
ucts of each SUSY particle contain an invisible y?, no mass peaks can be reconstructed 
directly. Instead, masses must be inferred P| from kinematic endpoints and other prop¬ 
erties of the events. Developing methods to do this has been a main emphasis of the 
studies to date. Typically, events are simulated for a particular SUSY model and for 
the Standard Model backgrounds using a parton shower Monte Carlo program such as 
HERWIG [|], ISA JET 0, or PYTHIA 0, the detector response is simulated using a fast 
parameterization, cuts are made to give a good signal/background, and various kinematic 
distributions are reconstructed. A number of examples are presented below. 

2 Search for SUSY 

Since g and q are strongly produced, their cross sections are comparable to QCD at the 
same Q^. UR parity is conserved, their decays produce distinctive events with large 
IpT- A typical analysis requires at least four jets with Et > 100,50,50,50 GeV and 
IpT >100 GeV and plots as a measure of the quantity 

Tfeff = IpT + Y PTj ■ 
jets j 

For large Meg the Standard Model background is typically 10% of the signal. 



mo (GeV) 


Catania 18 


Figure 1: Plot of 5cr reach in jets + I^t channel for mSUGRA model |^. 


The reach is thus limited mainly by the g and q production cross sections and hence 
masses. This reach is shown for the minimal SUGRA model in Fig. 2. Note that the 
sensitivity for one month at 10% of design luminosity is about 1500 GeV. Of course it 
will take more than one month to understand the detectors, but nevertheless one might 
hope for an interesting talk at SUSY08. Squark and gluino decays typically involve one 
or more intermediate charginos or neutralinos, giving rise to multi-lepton signatures as 
well as ones involving jets. The reach for such signatures is shown in Figure]^. 

The emergence of the SUSY signal from the Standard Model background is illustrated 
for a typical case in Figure ^ Also shown is the correlation between the peak of the Mgs 
distribution and 


AfsuSY = 


Ei Mjaj 




SUSY 


MHx^) 

AfsusY 


Obviously, however, any such inclusive correlation is model dependent. 

While the reach plots in Figure apply only to mSUGRA, similar reach in gluino 
and squark masses should apply to any model in which they decay into an invisible and 
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Figure 2: Reach limits in various channels for 100 fb ^ |^. 



Figure 3: Left: Typical distribution. Right: Correlation of Meg with SUSY mass 
scale . 















Figure 4: Plot of e’''e , fi^fi and mass distributions for mSUGRA with 2-body 

cascade decay (left) and Z plus 3-body decay (right) 

relatively light Xv Some models are easier. For example, GMSB models with prompt 
Xi —>■ Gx or i ^ Gi decays provide additional photon or lepton handles to suppress 
Standard Model backgrounds. GMSB models with a quasi-stable i provide a penetrating 
charged particle with P < 1 that is well measured by the muon detectors. AMSB models 
give few single leptons but larger ]pT- If R parity is violated via y? —>• t^l~v or qql^ qqp, 
there are again additional leptons. Perhaps the most difficult case is i?-violation via 

—>• cds, giving multiple jets but no b jets. The background is poorly known, but it 
seems likely that one must rely on leptonic cascade decays. In general, however, discovery 
of gluinos and squarks with masses of order 1 TeV appears straightforward. Masses of 
3 TeV are of course much more difficult and would probably need a luminosity upgrade. 

3 SUSY Measurements using Leptons 

If R parity is conserved, then the Xi is invisible and there are no mass peaks. Nevertheless, 
kinematic endpoints allow one to determine mass combinations. The simplest example is 
the decay 

For the dilepton signature, events are selected to have two leptons with (typically) 
Pt/ > 10 GeV and \rii\ < 2.5 in addition to multiple hard jets and large IpT- Then 
the dominant Standard Model background is tt. The background from ti or from any 
other Standard Model source involving two independent decays cancels in the combination 
e+e“ -|- The same is true for two independent chargino decays in mSUGRA, 

and it is likely to be true in any SUSY model that avoids p ^ ex and other low-energy 
constraints. 

An opposite-sign, same-flavor dilepton signature can arise from a 2-body cascade decay 
X 2 —^ —*■ from a 3-body x^ decay, or from a decay via a Z. All 

three cases are illustrated in Figure They are clearly distinguished by shape. Such 
endpoints are observable over a signihcant fraction of the mSUGRA parameter space, as 
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Figure 5: Reach for observing dilepton endpoints in SUGRA models with 1 fb lOfb ^ 
and 100 fb“^. Theory (TH) and experimental constraints are also indicated [Q. 



Figure 6: Dilepton + jet distributions for mSUGRA Point 5 as described in the text. 


illustrated in Figure |^. In particular, a large part of the mSUGRA parameter space with 
acceptable cold dark matter has light sleptons and hence enhanced decays. 
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Figure 7: Left: Scatter plot of reconstructed values of mi = vs. mi = M^o for LHC 
Point 5 (S5) and for an “optimized string model” (01) using multiple measurements from 

Projection onto M^o axis []H 


the decay chain —>■ X 2 Q Q- The stars mark the input values. Right: 


When a longer decay chain can be identified, more combinations of masses can be 
measured. For example, at mSUGRA Point 5, the dominant source of dileptons is qi 
X 2 Q Xi^^^~Q- The gluino is heavier than the squarks, and the hardest jets 

generally come from the squarks. In such a case, one can form an i~^£~q endpoint and 
two £^q endpoints in addition to the £^£~ one P, In addition, if a lower limit on the 
£^£~ mass is imposed, there is also an £^£~q threshold. 

All of these distributions after experimental selections are shown in Figure p. The 
endpoints and thresholds can all be expressed in terms of the masses involved using 
elementary kinematics, and they provide enough constraints to determine the masses. 
The results taking estimated errors into account for the Xi f/j masses are shown in 
Figure!^ for two models with similar masses. All the masses involved are determined quite 
accurately as functions of M^o; in particular the two models can be easily distinguished. 
The Xi mass is determined to about 10% through its effect on the kinematics even though 
it is invisible. 

An alternative approach is to reconstruct X 2 momentum assuming that the Xi mass 
is known and using 


Px°2 = Pee 


1 + 


Mfp 


While this is exact only for a 3-body X 2 fhe endpoint, it can be used as 

an approximation elsewhere. An analysis for SUGRA Point B (mo = 100 GeV, mi /2 = 
250 GeV, Ao = 0, tan /3 = 10, > 0) using this method looks for g ^bb ^ ^ ££bb 

Xi££bb. The dilepton distribution after cuts and the reconstructed x^b distribution using 
75 < Mil < 92 GeV are shown in Figure A second b can be added to reconstruct the 
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Figure 8: Left: Dilepton distribution for Point B analysis. Right: M{)(^h) distribution 
using inferred X 2 niomentum||15||. 


gluino mass; see Ref. 
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4 SUSY Measurements using h —> hh 

If the decay X 2 is kinematically allowed, it often has a substantial branching ratio. 

The decay h ^ bb can be reconstructed by selecting events with multiple jets plus large 
]pT and plotting the mass of pairs of jets tagged as 6’s. The reach for this signature in 
mSUGRA, shown in Figure |^, covers a substantial fraction of parameter space. It might 
well be the discovery mode for the light SUSY Higgs. 

If /i —> 66 can be reconstructed, it can then be combined with jets to determine other 
masses. In Figure |I^ for example, the h was combined with each of the two hardest jets 
in the event; the smaller of these masses should be less than the q ^ ^hq endpoint. 
This measurement is less precise than those involving leptons, but it may be useful if the 
leptonic branching ratios are small. 


5 Complex SUSY Signatures 


SUSY events can be much more complex than those for the cases discussed above. An 
example is the “focus point” region: in mSUGRA as mo becomes large for hxed values 
of the other parameters, —>■ 0 and radiative electroweak symmetry breaking fails. Near 
this boundary, the dominantly Higgsino, leading to acceptable values of cold dark 
matter [^. The exact location of the focus point region is very sensitive to details of the 


SUSY spectrum calculation |T^, but such a region surely exists 

m 


In one calculation 


the mSUGRA point mo = 1500 GeV, mi /2 = 300 GeV, Aq = 0, 
tan/3 = 10, /i > 0 is close to the // = 0 focus point boundary. Since the squarks are 












































tan (3 = 2, n < 0 tan (3 = 2, ^ > 0 




tan |3 = 10, n < 0 


tan p = 30, n < 0 




Figure 9: Plot of bb dijet mass distribution (points) with h ^ bb signal (solid), SUSY 
background (dashed), and Standard Model background (shaded) for various tan/3^^ 


heavy, SUSY production is dominated by gg pairs. The largest g branching ratios are 
calculated to be B{g —>■ Xitb + h.c.) ~ B{g —>■ y ;2 + h.c.) ~ 23%. These decays give 

events with 12-16 jets and leptons! The methods described previously suffice to find such 
signals. Many of the basic starting points discussed previously also work; see for example 
Figure 0. But sorting out the combinatorial background in such complex events is quite 
difficult and has not yet been solved. 

i?-parity violation with Xi QQQ, and especially Xi cds is another example of a 
SUSY model with very complex signatures. In such models the signal has very high jet 
multiplicity — nominally six jets just from the y? decays — few or no b jets that might 
be used to reduce the combinatorial background, and small IpT- The QCD background 
for such events is not well known; the 1-loop QCD corrections to the production of 6-10 
jets seem unlikely to be calculated even by the end of LHC running. Estimates based on 
QCD shower Monte Carlo programs suggest that the backgrounds are comparable to the 
signals, but this is at best a rough estimate. 








Figure 10: Minimum M{hj) mass from h ^ bb combined with two hardest jets 0. 
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Figure 11: Dilepton mass spectrum for focus point SUSY model[^ pT| . 


e.g. 


For such scenarios, therefore, it is probably necessary to rely on SUSY cascade decays, 

qqqi^i ~, to provide additional handles to suppress the QCD 


X2 


background. One example that has been studied uses such decays with branching ratios 
taken from mSUGRA Point 5. Even though there are no missing particles, reconstruction 
is difficult because the jets from the y? soft and the mass resolution is only ~ 10%. 
Dehne the mass combinations = M{qqqi^i~) ± M{qqq)] then the jet resolution 
largely cancels for M~, as can be seen in Figure |^. The distributions in Figure IR show 
reasonable distributions in these variables. 

These are just two examples of SUSY scenarios giving complex signatures. Such sce¬ 
narios are certainly possible, and much more work is needed to develop analysis strategies 
to deal with them. 
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Figure 12: Scatter plot of M+ vs. M for violating x? —^ QQQ decays [p2 



Figure 13: 
events with 


Distributions of qqq) ± M{qqq) in 

Xi QQQ after selecting either the peak or sidebands in 


i?-parity violating SUSY 

M- 


6 SUSY Signatures with r’s 

The initial motivation for studying r signatures came from the fact that for large regions of 
mSUGRA parameter space with tan/? S> 1 the only two-body decay modes are X 2 
and xf Then these modes have branching ratios close to 100%, and all SSUY 

decays involve r’s. Such scenarios are disfavored because they tend to give contributions 
to < 7 ^ — 2 much larger than observed. 

Low-energy tests of the Standard Model seem to require that the hrst two generations 
of squarks and sleptons be nearly degenerate. This degeneracy is not understood in 
general but does occur in mSUGRA. Even in mSUGRA, however, the third generation of 
squarks and sleptons is split by renormalization group and left-right mixing effects. Thus, 
study of the third generation is likely to be essential for understanding SUSY. 

Because of technical constraints, the ATLAS and GMS vertex detectors do not have 
sufficient resolution to identify r decays cleanly. Hence, leptonic r decays cannot be 
distinguished from prompt leptons, and r’s can only be identihed using narrow, one-prong 

































Figure 14: QCD jet rejection vs. efficiency for hadronic r decays based on full simulation 
of the ATLAS detector |^. 



Figure 15: Plot of visible mass distribution with hadronic r decays at LHC 

SUGRA Point 6 §. 


jets, for which there is substantial QCD background. As can be seen from Figure the 
typical jet rejection is a factor of about 100 for a r efficiency of 50%. 

For mSUGRA with mo = mi /2 = 200 GeV, Aq = 0, tan/? = 45, the only 2-body 
decays of light gauginos are y® ^ind xf These decays therefore dominate. 

Since there is ]pT both from z^^-’s and from x?’s, the r distributions can only be inferred 
from their hadronic decay products. Instead of a sharp rr edge at 59.64 GeV, one would 


observe the distribution shown in Figure 15 


The visible r momentum depends both on the momentum and on the polarization of 
the parent r. The decay r —is maximally sensitive to the polarization while high- 
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Figure 16: Left: Reconstructed masses for various r decay modes. Right: rr mass 
distribution for all modes, for one r —>■ vrz/ decay, and for decays with reversed helicities. 


mass decays are rather insensitive to it. A study of the separation of different decay modes 
has been carried out using a full GEANT-based simulation of r’s from SUSY events in 
the ATLAS detector. The visible r mass is reconstructed by combining charged tracks 
and electromagnetic calorimeter cells, exploiting the hne granularity of the calorimeter. 
The reconstructed masses. Figure |T^, show a good separation of nv decays and some 
difference between pv and aif decays. The same hgure indicates the sensitivity of the rr 
mass distribution to a reversal of helicities. 


7 GMSB Signatures 

In GMSB models the lightest SUSY particle is the light gravitino G. The phenomenology 
depends on nature and lifetime of the next lightest particle (NLSP), which can be either 
the Xi or S' Decays into G lead to more distinctive events and to longer decay chains, 
such as X 2 ^ Gxi^i~, that can be used to determine SUSY masses |^. 

The lifetime for the NLSP to decay to a G is unknown and is an important parameter 
of the model. Short Xi Gj lifetimes can be measured by selecting Dalitz decays and 
using the vertex detector. For cr^o 3 > 1 m the signature is occasional photons that do 
no point to the primary vertex. The ATLAS electromagnetic calorimeter was designed 
to measure photon angles with A6 ^ improve the h ^ 77 mass resolution. It 

also has good time resolution. At ~ 100 ps. Both can be exploited to detect long-lived 
Xi —>■ Gj decays with good efficiency, as seen in Figurefeftdr20-65. If no signal is seen, 
an upper limit of order cr ~ 100 km could be established |^. 
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Figure 17: Angle (left) and time delay (right) for photons from long-lived Xi 
the ATLAS detectorQ. 



Figure 18: Left: Time of flight separation of muons and quasi-stable i Right: Mea¬ 


surement of £ lifetime for different estimates of the systematic error on the acceptance 


A long-lived slepton NLSP would give tracks with (3 < 1 through the calorimeter and 
muon system. These can be well identihed by using the muon system as a time of flight 
system, as can be seen in Figure |T^ p, The lifetime can be estimated by counting 


the number of events with zero, one, and two such tracks, as is also shown in Figure |I8 


It should also be possible to determine the lifetime by looking for kinks in the central 
tracker, but developing the needed pattern recognition algorithm is non-trivial and has 
not yet been attempted. 





















8 Outlook 


If SUSY exists with masses below about 1 TeV, it seems likely that the ATLAS and 
CMS detectors should find evidence for it shortly after the LHC begins operation. Only 
a limited number of SUSY models and cases have been investigated so far, and in all 
cases the answer has been known. These studies seem sufficient, however, to sketch the 
broad outlines of an initial program of SUSY measurements (assuming that R parity is 
conserved): 

1. Search for an excess of multijet + JpT events over Standard Model predictions, which 
of course must be checked against other measurements. 

2. If such an excess is found, select a sample of SUSY with simple cuts such as those 
described above. 

3. Look for special features such as y’s or long-lived in these events. 

4. Look for h jets, hadronic r decays, etc. 

5. Try simple endpoint-type analyses such as those described above. 

This program looks quite feasible. 

There is of course much more information available, including cross sections, branch¬ 
ing ratios, and other kinematic distributions. Ultimately one will want to use the hrst 
measurements to guide further analyses incorporating all the information. This program 
will need the full power of the LHC and its detectors. 
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